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The behavior of a long Josephson tunnel junction drastically depends on the distribution of the 
dc bias current. We investigate the case in which the bias current is fed in the central point of a one- 
dimensional junction. Such junction configuration has been recently used to detect the persistent 
currents circulating in a superconducting loop. Analytical and numerical results indicate that the 
presence of fractional vortices leads to remarkable differences from the conventional case of uniformly 
distributed dc bias current. The theoretical findings are supported by detailed measurements on a 
number of 5-biased samples having different electrical and geometrical parameters. 

PACS numbers: 03.70.+k, 05.70.Fh, 03.65.Yz 



Symmetry principles and how they are broken are fundamental concepts in physics. A recent experiment ^ |2] 
demonstrated that the symmetry can spontaneously break during the fast superconducting phase transition of a metal 
ring and both fluxoids or antifluxoids can be trapped in the ring while its temperature crosses the superconducting 
critical temperature. This phenomenon was predicted a long time ago as one among several possible condensed matter 
cosmological experiments [3, to check the validity of the causality principle in the early Universe jl]. In the experiment 
of Ref.[JJ the presence of the persistent currents associated with the flux trapped in the ring relied on the radial 
magnetic field modulation of the critical current of a planar Josephson tunnel junction (JTJ) having the peculiar 
configuration shown in Fig. [TJa). It consists of a ring shaped junction cut at some point with the bias current fed at 
the diametrally opposite point. The cut leaves a gapped annular junction and reliefs the junction from the constraint 
of the 27r-periodic boundary condition of annular junctions 5 . In that experiment the ring itself acted as the junction 
base electrode, while the top electrode had the shape of a circular arc of about 300°. Later on it will be demonstrated 
that a gapped annular junction in a radially uniform magnetic field H r is topologically equivalent to a linear junction 
in an in-plane uniform field Hu as that depicted in Fig. fl[b). 




The task of this work is to study the properties of a S-biased or single point injected overlap Josephson tunnel junction 
whose physical length is L, i.e., — L/2 < X < L/2 and whose width is W, i.e., — W/2 < Y < W/2. To simplify the 
analysis, we assume that Josephson current density J c is uniform over the barrier area and that the junction width W 
is smaller than the Josephson penetration depth Xj. Ideally, the bias current is fed to the junction by infinitely narrow 
electrodes. In real devices the <5-bias approximation is achieved as far as the electrodes carrying the bias current in 
and out of the tunnel barrier are much narrower than the junction Josephson penetration depth Xj] however, for 
very long junctions, it is only required that the electrodes widths are much smaller than the junction length L. In 
window type JTJs one more requirement is that the passive region surrounding the tunnel area, the so-called idle 
region [HI [5], needs to be narrower than the current-carrying electrodes otherwise the bias current diffuses before 
entering the barrier and the sharp bias profile gets smeared. We will consider both intermediate length (L ~ Xj) and 
long (L >> Xj) JTJs (the behavior of small junctions is not affected by the bias profile). We remark that with the 
current injected at the junction extremities we recover the well-known case of so called in-line configuration treated 
by pioneering works on long JTJs soon after the discovery of the Josephson effect [HI EH E] ■ For in-line JTJs it is 
important to distinguish between the symmetric [TUJ E] and asymmetric configuration ITU]: the former is achieved 
when the bias current enters at one extremity and exits at the opposite one, while the latter is obtained when the 
bias current enters and exits from the same extremity. Since in this paper we will only consider the case in which 
the bias current is fed in the middle of the junction long dimension, for symmetry reasons, we do not need to specify 
the electrode configuration. Nevertheless the more general situation in which the bias current enters and leaves in 
two generic lateral points along the junction will deserve consideration in the future. We will show that, with the 
bias current centrally injected, both the static and dynamic junction properties reveal interesting phenomena whose 
understanding will serve as a base for the study of more general cases. In this regard, a considerable attention has 
been recently given to the case of — tt transition Josephson tunnel junction obtained by closely situated current 
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FIG. 1: (a) gapped annular and (b) linear 5-biased Josephson tunnel junctions. The base electrodes are in black, the top 
electrodes are in gray and the junction areas are white. A gapped annular junction in a radially uniform magnetic field H r is 
topologically equivalent to a linear junction in an in-plane uniform field Hu. 

injectors[TU [13]. For the sake of completeness it is worth to remark that the lateral point injected bias was already 
treated in the literature, but only in the limit of infinitely long JTJs |14l IT51 IT5] . 

THE MODEL 

It can be shown that for a (5-biased JTJ, as that sketched in Fig. [2] the gauge-invariant phase difference <j> of 
the order parameters of the superconductors on each side of the tunnel barrier obeys the static or d.c. perturbed 
sine-Gordon equation: 



A 2 d 2 4> 
j dX 2 



in which the term 7 = h/Io is the external bias current lb normalized to Iq — J C LW and 6 is the 5- function, 
normalized units of x = X/Xj, the above partial differential equation (PDE) becomes: 



In 



dx 2 



= sin < 



(x) — j8(x). 



(1) 



If the time t replaces x in Eq.Q and a new phase tp = 4>— n is defined, we obtain the equation for a lossless pendulum 
with an applied torque like a train of pulses with unitary period. In this framework, the boundary conditions will 
become constraints on the pendulum angular velocity half period before and after each pulse. 

It is possible to derive that, introducing the Heaviside step function TL{x), the Lagrangian density C of our system is 
constant: 



Y+cos</> + 7<fe(0)«(z). 
By subtracting the constant quantity j<f> x (0)/2, we can rewrite the last expression as: 



| + cos0 + ^#sgn(,) = C, 



(2) 
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FIG. 2: 3D sketch of a centrally injected <5-biased linear planar Josephson tunnel junction. 



with C being a constant depending on both the external bias and magnetic field. Note that if <p x (x) is discontinuous 
in x = 0, then 2<f) x (0) = X (O+) + cf> x (0— ). The first two terms in Eq.Q are related to the system free energy density 
£(x) = 1 — cos0 — <fi 2 ,/2, while the last term represents the two-level potential U(x) generated by the 5-shaped forcing 
term. 



Boundary conditions 



The magnetic Josephson equation[T7] states that the phase gradient is proportional to the magnetic field: 



V</>=— Hxn, (3) 

<3>0 

where n is a unit vector normal to the insulating barrier separating the two superconducting electrodes. If the two 
superconducting films have thicknesses and London penetration depths \li,2 and tj is the barrier thickness, then 
the effective magnetic penetration d e is given by |18j : 



d e = tj + Aii tanh — h X L2 tanh — — , 

which, in the case of thick superconducting films (di » Xu), reduces to d e w Xli + Xl 2 (since always >> tj). 
From Eq.(|3|) it follows that, for a linear JTJ in an external uniform field Hu applied in the junction plane perpendicular 
to its length L, the boundary conditions are: 



dX 



X=±L/2 



(4) 



with k = 2nd e fio/$Q. Being X~ 2 — 2Trd e noJ c /&o — K Jc (3>o is the magnetic flux quantum and fj,o is the vacuum 
permeability), we have nXj = l/J c Xj. Introducing the critical field H* for a unitary junction with a Fraunhofer 
magnetic diffraction pattern H* — (fiod e Xj) = 2irJ c Xj = 2tt / nXj , we get nXj = 2tt/H*. In normalized units of 
h = nXjH\\ = 2ttH\\/H* 7 the boundary conditions (E| for Eq.dTl are: 



dx 



x=±l/2 



(5) 
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in which we have introduced the junction normalized length I = L/Xj. Note that, with this notations, the normalized 
critical field h* of a short JTJ equals 2?r/l. For the gapped annular junction in a radially uniform field H r , the 
boundary conditions for the PDE in Eq.([T|) are independent on the gap angle [19] and coincide with those in Eq.pJ, 
but now with h = 2ix H r /H*. In other words, a linear junction in an in-plane field H\\ and a gapped annular junction 
in a radial field H r arc governed by the same PDE with the same boundary conditions. Therefore, for the remaining 
of the paper we will use the properly normalized field h both for linear and gapped annular JTJs. A radial magnetic 
field can be generated by a current flowing in a control line in the shape of a loop concentric to the annulus; however, 
the simplest wayp] is to have a ring shaped base electrode and to apply an external field perpendicular to it to 
induce tangential screening currents proportional to the applied field. The only disadvantage of this method is that 
the effective radial field felt by the gapped annular junction depends on geometrical factors such as the ring inner 
and outer radii and the junction position relative to the ring. The radial dependence of the current circulating in a 
superconducting ring has been calculated under many different conditions [2"U] . 

From Eqs.pj) and (|5|, we have that the junction Josephson current ij is proportional to the bias current; in fact, 
being J" b[x)dx = 1, we get: 



i ^| = l/_ i / /2 sin ^ )cfa= T(S 

that is, for a fixed bias current the zero- voltage current Ij passing through a <5-biased JTJ is inversely proportional 
to its normalized length (Ij — I^Xj/L). In the well-known case of an overlap JTJ with uniform bias 7(2;) = 7 tt , it 
would be ij — 7 U (i.e., Ij = It), meaning that 7„ cannot exceed unity. In contrast, Eq.([6]) implies that for (5-biased 
junctions the largest value 7 C that the normalized bias current can achieve is determined by the junction normalized 
length. Later on it will be found that, in our case, I7I < 4. 

The jump in the phase gradient at the injection point (the axis origin in case) can be calculated directly from Eq.([T]). 
For any < x < 1/2 we can write: 



M/2 



dx 



x=-l/2 1 



1 = 1 



(6) 



dx 



dx 



[sin 4>{x) — r yS(x)]dx = / sin (f>(x)dx — 7 



Taking the limit x$ — > 0, the integral vanishes, being 4>(x) a continuous function; then the phase derivative jumps at 
the injection point x — 0: 



d4 

dx 



x=0+ 



dx 



(7) 



This equation was first reported by Kuprianov et al. 
injection. We can rewrite Eq.Q as: 



[see also Ref. [T5] at (par. 8. 5)] for JTJs with lateral single-point 



dx 



dx 



7/2- 



(8) 



x=0+ 



in which the constant ho is a measure of the (exponential) penetration of the external magnetic field (if any) into the 
center of the junction and, in general, cannot be determined a priori. 



ANALYTICAL APPROXIMATE SOLUTIONS FOR LONG JTJS 



In this paragraph we will discuss the possible analytical approaches to determine the phase profile of a point 
injected long JTJ in the Meissner state. In absence of external bias and magnetic field the Josephson phase profile 
4>(x) is identically equal to zero. If 7 and h are sufficiently small, then |0(x)| << 1, so that Eq.(JT|) can be linearized, 
being smcj) Ri cj). Barone et al. [5TJ pointed out that a piecewise linear current phase relationship jj = x 2 ^ < 1 
(with x G [0-5,1]) can correctly handle the case when |</>(x)| < n/2. Obviously in a linear approximation the phase 
profile looses its correct structure within the scale of Xj, but averaged results such as the Josephson current are not 



5 



3.5 
3.0 
2.5 
2.0- 
I 5 
I 
0.5 
0.0 



\>20, h=0, y=4 




2.0 

1.5- 

1.0 

0.5 

e 

X o.o 

-0.5 
-1.0 
-15 
-2.0 



-10 -8 



(a) 



-2 
X 

(b) 



8 10 





FIG. 3: (a) Phase profile <j) e {x) for an infinitely long 5-biased JTJ as in (HO) with £ = 0, for x G [-10, 10]. (b), (c) and (d) 
show, respectively, (f>%{x), 4>xx( x ), an d sin0 e (:r). It is cos(f> e (x) = 1 — [4>%(x)\ z /2. 



significantly affected by the approximation. However, as will be shown later on, for <5-biased JTJs in the Meissner state 
it is |0(x)| < 7r so that also the piecewise linear approximation fails and a cubic approximation can be conveniently 
adopted sin<^> ~ a<f> — ftcj) 3 leading to a Duffing-like differential equation for the oscillation of a soft (a, (3 > 0) spring 
system. In the special case of the underdamped and unforced Duffing equation, exact solutions can be written in terms 
of Jacobi's elliptic functions [23 . Further, approximate solutions of the forced Duffing equation could be found using 
the perturbation methods [2JJ when f3 « a. However, simple solutions can be readily obtained in the approximation 
of very long junctions. 

h = and 7 / 0, 1/2 > 2n 

We begin with the case of no externally applied magnetic field which, in the well-studied case of uniform bias, results 
in a constant phase profile 4>{x) — sin -1 7„ (mod 2tt). For symmetry reasons, in the case of 5-biased JTJs, the solution 
of Eq.Q has to be an even function 4> e {x). Therefore its gradient is an odd function such that 4>%{ x ) + 4>%{— x ) = 
for any x e]0, Z/2], meaning that ho = in Eq.Q. With ho = 0, Eq.Q provides two extra conditions on the phase 
left and right derivative at the injection point: 



dx 



dx 



7 /2. 



(9) 



x=0+ 



For I >> Ait (strictly in the limit I — > 00) with the above conditions, an approximate solution of Eq.Q is a cusp- like 
function: 



<f> e (x) = ±4tan" 1 exp-(|a;| + £), 



(10) 
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in which £ is a non-negative constant set by the bias current 7; Eq.([To| can also be cast in the form[5]: sin </> e (x)/2 = 
±sech(|a;| + £). It is cf) e (±oo) = </>|(±oo) = and (f) e (0) = ±4tan _ ^~exp(— ^). The sign in front of Eq.([To| concords 
with that of the external bias current. With £ = (and positive 7), then 4> e (0) = tt, meaning that Eq.( 10 1 corresponds 



to a scmifluxon (tt jump) on the left side of the junction and anti semifluxon (— tt jump) on the right side, as shown in 
Figs. [3^a-d) for —10 < x < 10 (we observe that, as required, 4>xx — sin^>, everywhere, but for x — 0). The dependence 
of £ on 7 can be found observing that, in force of Eq.Q, 7 = (j>%(0— ) — <j)%(0+) = ±4sech£, so the largest possible 
amplitude value for the normalized bias |7 C | — 4 is achieved when £ = 0. For £ < the phase in the origin grows above 
the threshold value |0 e (O)| = tt and the static solution in Eq.([To| is no longer stable: the two semifluxons develop 
into integer fluxons driven away in opposite directions under the effect of the Lorentz force associated with the bias 
current. According to Ref.[15 and from numerical simulations (reported later), we found that the instability arises 
when the amplitude of 7 exceeds the critical value j c = 4, corresponding to Ic = ±4J c XjW. With such a notation 
7 = ±7 c sech£, i.e., £ = cosh _1 (7 c /|7|). The last expression allows us to find the dependence of 4> e (0) on 7; it is 
found that cos0 e (O) = 1 — 2(-f/j c ) 2 (i.e., sin0 e (O)/2 = j/"f c )- We remark that, inserting the value 7 C = 4 in Eq.Q, 
we reach the important conclusion that for very long junctions the normalized zero-field critical current is inversely 
proportional to the the junction length: 



i c (h = 0) 



I 



(11) 



as for asymmetric inline junctions^ with the only difference that j c = 2. Indeed, a long <5-biased JTJ in zero 
external field is equivalent to two inline asymmetric junctions in a parallel configuration; this will not any longer 
be true in presence of an external field. For a generic 7 value in the interval [—4,4] the phase difference A(j) 
'/;:(). — 00) = (f)(0) — 2 sin -1 7/7 c corresponds to what in nowadays language is called a /c-fractional vortex where 
k — A0/(2-7r) = 7T _1 sin -1 7/7 c . Semi (£ = 0) and fractional (£ > 0) vortices are presently receiving a great deal of 
attention in the context of — tt transition Josephson junctions[12 25 26 27 . 

In summary, the main difference between short and long <5-biased JTJs is that in the former case the solution becomes 
unstable when 7 > j c — I and 4> e (0) > tt/2, while in the latter case the phase profile becomes unstable when 7 > j c = 4 
and 4> e (Q) > tt. The gradual crossover from intermediate to long JTJs has to be calculated numerically; it was found 
to be nicely described by the following empirically found relationships: 



7 C (ft = 0,1) = 4tanh 



/ 

4' 



(12) 



(j> e c (x = 0, ft = 0, 1) = — + tan 1 cxp(^ - tt). 



(13) 



hj^O and 7 = 0, 1/2 > 2n 
With 7 = 0, the nonlinear PDE Eq.Q reduces to: 



dx 2 



sin 4>(x). 



(14) 



Eq.(14l was first introduced in the analysis of long asymmetric inline JTJs by Ferrel and Prangc[9 in 1963; however, 
four years later Owen and ScalapinofIT] reported an extensive study of its solutions for long symmetric inline JTJs: 
this is why Eq.(14l is commonly known as the equation of Owen-Scalapino(OS). With boundary conditions as in 
Eq.(|5]), the solution of the OS equation has to be an odd function <p°(x) [i.e, with 4>° x (x) even]. For given h < 2 and 
/ > 7r/2, exact solutions exists in terms of Jacobian elliptic functions. However, upon the assumption that the JTJ is 
so long that the magnetic field in its center can be neglected, a simple approximate solution exists. In fact, for / >> 1 
(in practice I > Air), the solution of Eq.(14l with 4>°(0) = is: 



b°(x) = ±4 



tan 1 exp [ x 



— tan exp —x + C + 



I 



(15) 
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FIG. 4: (a) Phase profile 4>°{x) as in Eq.p5| for 7 = 0, h -. 
(d) its cosines. Note that sin0°(:r) = (f>^. x Jx) an d cos(j> (x) 



2 and I = 20; (b) its first derivative, (c) its second derivative, and 



in which the non-negative constant C, is set by h = </>°(±Z/2) w ±2sech£ indicating that the largest possible amplitude 
of the normalized magnetic field[15] is h c = 2 corresponding to H c = 2J c Xj and £ = 0. Now <fi o (0) — ±4[tan _1 exp(C + 
1/2) — tan -1 exp(C + 1/2)] — 0. The sign in front of Eq.( 15 1 now concords with that of the applied field. Eq.( 15 1 is 
shown in Figs. |4]for I = 20 and h — h c . Again Eq.(15l, as Eq.(10), represents a superposition of two static fractional 
vortices, but now they a pinned at the junctions extremities and not in its middle point. As \h\ exceeds h c , i.e. £ < 0, 
we exit the Meissner regime and some magnetic flux enters into the junction interior; one (or more) integer vortices 
(fluxons) gradually develops at each extremity and move towards the center resulting in a phase profile that can no 
more be written in terms of Jacobian elliptic function. For h » h c the phase profile resulting from the superposition 
of several closely packed fluxons will be approximately linear <j)°(x) — hx and we recover the behavior of small JTJs 
(Fraunhofer regime). 

The gradual crossover from intermediate to long JTJs has to be calculated numerically; here we anticipate that it was 
found to be nicely described by the following empirical relationship: 



2tt I 
h c {l) = -r + 2tanh— , 

in which the first (Fraunhofer) term dominates for small values of I, while the second (saturating) term dominates 
for large I values. It is worth to remark that if no current feeds a junction, then its electrode configuration does not 



affect the phase profile; in other words, inline, overlap and ^-biased JTJs all have the same phase profile as in Eq.(15l 
in presence of a given external magnetic field h. 
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h 0, 7 0, 1/4 > 2tt 

Approximate static phase profiles for non-zero h and 7 values are obtained observing that Eq.Q can be rewritten as 
two identical, although independent, OS PDEs for the left (—1/2 < x < 0) and right (0 < x < 1/2) inline asymmetric 
half-junctions each with an effective bias current 7/2. The boundary conditions for the left and right half junctions 
are simply given by Eqs.([5]) and (pj: 



dx 



-1/2 



dx 



dx 



dx 



=1/2 



h. 



(16) 



(17) 



It is convenient to rewrite the constant C in Eq.([2| as C = (2 — k 2 )/k 2 \ for C > 1, then k 2 < 1 and vice-versa. 



According to the notation of Ref.[TT], the boundary conditions in Eqs.(|T6[) and (17 1 can be rewritten in terms of 
Jacobian elliptic functions dn(ir,fc^) and cn(x,l/k 2 ) of argument x and modulus k 2 d and 1/fc^, respectively. For 
k d < 1: 



while for k r > 1: 



(In 



1/2 - x d 



2 



and 



dn 



-X d 

kd 



2 2' 



cn(Z/2 



X c 



and 



cn (-x 0c , 1/kf) = ± 



2 2' 



Being the above mentioned elliptic functions limited to the [—1,1] range, the solutions of the OS problem can be found 
as far as both \h\ and \ j/2\ are smaller than 2. Once h € [—2, 2], 7 € [—4,4] and I > it are given, the couples (xod, kd) 
or (xoc, k c ) can be numerically found. This mathematical procedure allows to find many (sometime physically non- 
interesting) solutions. It is well known that the number of possible solutions increases with the junction normalized 
length [25]. 



For 1/2 » Air, the approximate solution of Eq.(14) is: 



= ±4 
= ±4 



tan 1 exp (x — £) — tan 1 exp y— x + ( — - 
tan -1 exp {—x — £) + tan -1 exp [ x — ( — - 



for -l/2<x<0, 
for < x < 1/2 



(18) 



in which £ and £ are two non-negative independent constants; in fact, 7 « ±4sech£ and h w ±2sech£. If ft, = 0, then 
C — * 00 and 4>(x) — 4> e (x). Vice-versa, if 7 = 0, then £ — > 00 and </>(a;) = <j)°(x). In other words, the generic static 
phase profile of a long 5-biased junction in the Meissner state is obtained simply by the sum of (four) non-interacting 
fractional vortices. The sign in front of Eq.(18l concords with that of the product hj. Looking at the Eqs.(16) and 
(17 1, it is seen that an interesting situation occurs when h = ±7/2. For h = the bias current 7 flows symmetrically 
in the left and right junction sides [see Fig. |3jd)]. With h =/= the symmetry is broken and the current flows mainly 
in one of the junction sides. When \h\ = 7/2 < 2, then £ = £ and the applied current only flows in one junction side, 
while in the other side the average Josephson current vanishes. For h = 7/2 = 2 (£ = £ = 0) the expression above 
for <p(x) and its sine are depicted in Figs. [BJa) and (b), respectively. Only the junction right side contributes to the 
Josephson current. 

Let us remark that, being I7I < j c , in force of Eq.(|6]), for the largest Josephson current i c we have: 



(19) 



In other words, for long 5-biased JTJs in the Meissner state (—2 < h < 2) the critical current is independent on the 
externally applied field. There is no other junction configuration for which this peculiarity occurs. The result above 
is only apparently in contrast with Kuprianov[14 quadratic prediction for infinitely long point injected junctions: 



9 




h=2, y=4 



-10 -8 



6 8 10 



(a) 




FIG. 5: (a) <j>(x) as in ( |18| l for h = 7/2 = 2 (£ = £ = 0) and Z = 20; (b) sin0(a;). Only the junction right side contributes to 
the Josephson current. It is sia<j>(x) = 4> xx {x) and cos0 = 1 — [<^>(:r)] 2 /2. 



i c (h) = 



2 + \/4+7? 



(20) 



leading to i c (±h c ) = (1 + \/2)i c (0) /2 w 1.21i c (0). In fact, Eq.(20 ) also takes in the account the penetration of fluxons 



into the barrier leading to the non-Meissner regime characterized by higher modes phase profiles [TT] that will be 
discussed in the next Section. In 1985 Radparvar and Nordman|16j reported the experimental magnetic diffraction 
pattern of a very long (L ~ 100 Xj) Nb-Pb laterally injected JTJ in reasonable agreement with Eq.(p0| despite their 
electrode configuration only roughly realized the point injected approximation. 



STATIC NUMERICAL SIMULATIONS 



In this section we discuss the numerically obtained solutions of the PDE in Eq.Q with boundary conditions as in 
Eq.(|5]). The direct numerical integration of Eq.Q poses large problems of stability due to the fact that there are 
no losses in the system; to avoid this problems, we recurred to the integration of the time dependent perturbed 
sine-Gordon equation: 



4>xx - 4>tt -aia^)= j5(x) + afa - P<j> X xt, 



(21) 



with a = 3 in order to have a fast decay of the temporal features of the solution towards a static solution (in 
real device a < 0.01). The term containing the surface losses was simply dropped to save computer time, i.e. 
(3 = 0. Eq.(21l has been numerically integrated by using the commercial finite-element simulation package of Comsol 



Mutiphysics (www.comsol.com) for different values of the normalized length / which enters the PDE through the 
boundary conditions. The (5-function has been approximated by the continuous function fix) = psech 2 {2px) with the 
parameter p = O(10 2 ) [note that /(0) = p, /(1/p) = psech 2 2 « 0.07p and f(x)dx = 1]. 

In order to trace the different lobes of i c vs. h, it is crucial to start the numerical integration with a proper initial 
phase profile. As far as I < n the initial condition for the numerical integration of Eq.(JTJ was simply set to 4>{x) = 0. 
In order to find the several possible initial conditions for I > tt, the OS problem was solved for a given value of h 
and with 7 set to zero. Once the initial phase derivative was known, also the initial phase profile could be easily 
derived^ll . Finally, during the numeric calculation, 7 was changed until the numerical solution becomes unstable. 
We have numerically computed j c as the maximum allowed value of the bias current for each chosen value of the 
magnetic field h. Once "f c was found, the corresponding critical current i c could be calculated either numerically or, 
in a equivalent manner, resorting to Eq.(|6|. As expected, a general peculiarity of the numerically found phase profiles 
for long enough (5-biased junctions is that the phase values at the junction extremities {x — do not depend on 

7; viceversa, the phase behavior near the origin (x = 0) does not depend on the applied field h. 

Figs. [6^a)-(d) display the magnetic diffraction patterns for point injected JTJs with different normalized lengths. We 
remark that the i c {h) patterns are symmetric around h = 0. For comparison we also report the i c (h) for a uniform 
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FIG. 6: Numerically computed magnetic diffraction patterns for for different junction normalized lengths. From (a) to (d) i c 
vs. h respectively, for I = 2,4, 8, 16. The full dots refer to point injected bias current, while, for the sake of comparison, the 
open dots correspond to the well-known case of a uniformly biased junction. 



bias; more precisely, the full dots refer to the point injected current, while the open (gray) dots correspond to the 
uniform bias. For 2 = 1, the numerical data (not reported) closely follow the expected Fraunhofer-like dependence 
(with h c — 27r) with differences only in the third significant digit. As shown in Fig. [6]ja), for I = 2 we still have a 
Fraunhofer-like pattern. Pronounced deviations from small junction behavior were found for I = 4, but, as can be seen 
in Fig.[6jb), they disappear for large filed values. Increasing I, some ranges of magnetic field develop in correspondence 
of the pattern minima in which i c may assume two (or more) different values corresponding to different phase profiles 
inside the barrier. In fact, each pattern lobe is associated with a given vortex structure; more precisely, in the first 
lobe which, for I — 4, goes from h = to h c « 2.2, the external magnetic field is shielded and vortex cannot penetrate 
into the barrier (Meissner state). However, at the very end of this lobe a fluxon is present in the junction. In the 
successive lobes the magnetic field penetrates into the barrier and vortices are created in the barrier in a way which 
closely recalls the behavior of the type II superconductors, even though the vortices we are dealing with are quite 
different from the Abrikosov ones as they do not have a normal core. In the second lobe, moving from h ~ 3 to 
h w 4.6, we start from a phase configuration very similar to that at the right side of the first lobe in which one vortex 
is present in the barrier and we end up with two bunched fluxons. Adopting the terminology used in Ref. |llj. we refer 
to the first (Meissner) lobe as to '0 to 1 vortex mode' lobe, the second as the '1 to 2 vortex mode' lobe and so on. In 
general, one may talk about the 'n to n+ 1 vortex mode' when the junction contains more than n but less than n+ 1 
vortices. As I is increased from 4 to 8 drastic changes occur (the crossover point being approximately 2tt), as shown 
in Fig.^c): the critical current get smaller and smaller, according to Eq.(ll), and the principal lobe gets flatter and 
flatter resulting in a rather large plateau; further, the lobes broaden and overlap each other with h c converging to 2. 
This behavior is even more pronounced for I = 16 as shown in Fig. [6jcl) . In other words, the main (Meissner) lobe 
is a rectangle with corners in h — and 2 and 7 = and 4. The phase profiles in the non-trivial corners are those 
already shown in Figs. [3] [3]and[5j 

The phase profile corresponding to the 2 and 3 fluxon modes are reported in Figs. |7ja) and (b), respectively, for h 
(7) slightly below (above) its critical value. It is evident that the vortices only penetrate in left side of the junction. 
However, the situation reverses by reversing either h or 7. The maximum supercurrent associated to a given vortex 
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mode increases with the mode order; this is peculiar of 5-biased JTJs. 



EXPERIMENTS 



To investigated the properties of (5-biased JTJs, we have used are high quality Nb/Al — Al ox /Nb JTJs fabricated 
on silicon substrates using the trilayer technique in which the junction is realized in the window opened in a SiO 
insulator layer. We measured a large number of both linear and gapped annular junctions all of them with width 
W = 2 fj,m, but whose lengths varied from a L — 80 to 335 /im. The width of the electrodes carrying the current 
in and out of the barrier was 10 /im. The thickness of the S1O2 insulator layer was 200tito and the so called "idle 
region" , i.e. the overlapping of the wiring layer onto the base electrode was about 1 /^m for all the junctions. In order 
to vary the sample normalized lengths over a large range values, we used two sets of samples having quite different 
critical current densities (J c — 100 A/ cm? and 3k A/ cm 2 ), corresponding to Xj ~ 80 and ~ 12 ^m. In such a way 
samples were available with I spanning from about 1 to about 30. The values of the Josephson penetration depth 
were calculated taking into account the the effect of the lateral idle region[7l [8]. In Fig. [8] we report the log-log plot 
of the measured zero-field critical currents i c (0) of many 5-biased JTJs versus their reduced length L/Xj. The critical 
current has been normalized to the small junction theoretical value calculated as the 70% of the current jump at the 
junction gap voltage. We observe that, for large normalized length, the experimental data clearly follow an inverse 
proportionality law as expected from Eq.(lll and the transition from long to short junctions is nicely fitted by the 
empirical expression in Eq.( 12 1 (solid line). 

On real samples, the measurements of maximum supercurrent against the external field often yield the envelop of 
the lobes, i.e., the current distribution switches automatically to the mode which for a given field carries the largest 
supercurrent. Sometimes, for a given applied field, multiple solutions are observed on a statistical basis by sweeping 
many times on the junction current-voltage characteristic. Figs. [9ja)-(d) display the measured magnetic diffraction 
patterns for four linear samples in a uniform in-plane field with selected normalized lengths. Analogous results 
(not reported) were obtained for samples with the gapped annular geometry in an uniform, although uncalibrated, 
radial magnetic field. We found an excellent agreement with the results of numerical calculation discussed in the 
previous section. Marked deviation from the small junction Fraunhofer-like behavior were observed for L = 3Xj. 
It is evident that, for longer JTJs, the I C (H) is made of few flat segments corresponding to different static phase 
solutions containing more and more fluxons. Since the number of possible static solutions increases with the junction 
normalized length, it is straightforward to assume that, as I increases, the number of flat segments increases too, 
while their lengths decrease so that to merge in a monotonically increasing critical current as predicted in Ref. |14j 
and observed in Ref.[T^]. 

The dynamic of a (5-biased JTJ was also investigated both experimentally and by numerically integrating the time 
dependent Eq.(21l whose steady state solutions were typical of a soliton moving in a two-level potential, with the 
fluxon(s) being accelerated only when passing trough the origin. We have focused out attention on the zero-field single 



fluxon shuttling mode. Figs. 10 compares the first zero-field step profiles of a (5-biased (solid line)and a uniformly 
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biased (dashed line) long overlap JTJs having the same length (L = 150 /mi) and the same Josephson current density 
(J c = 3 kA/cm 2 ). The data were taken a T w 6iG Contrary to the case of uniform bias, a current threshold value 
is needed to avoid that fluxon stops due to friction. This value, of course, depends on the the junction length and 
losses: the longer the junction, the smaller should the losses be to maintain a finite-voltage dynamic state. The 
j-(4>t) plot, corresponding to the zero-field step profile, is smoother that that obtained when the bias is uniform and 
is characterized by fine structures associated to the fluxon interactions with plasma waves [121 130] originated by the 
fluxon itself when passing across the potential discontinuity. Often the back-switching transition has been observed 
at the top of the step, as indicated by the solid arrow pointing to the zero voltage state. As a matter of fact, point 
injected junctions seem not to reveal new dynamical states. We only remark that displaced linear slopes j25] were 
numerically found and experimentally regularly observed in absence of magnetic field and for junction with large 
losses (a w 1 or T > 0.7 T c ). when the d.c. bias exceeded its critical value. Displaced linear slopes were also observed 
in a presence of a magnetic field which, by increasing the field further, eventually develop in large amplitude flux flow 
steps. 



CONCLUSIONS 



In this paper we focused on <5-biascd (or point-injected) long Josephson tunnel junctions. Our analysis goes much 
beyond the previous theoretical [14 and experimental^^ works in which only extremely long junctions were treated. 
We considered a linear JTJ in a uniform magnetic field and a gapped annular junction in a radially uniform field. 
These two electrode configurations (shown in Figs, [ljb) and (a), respectively) are topologically equivalent, that is, 
they are described by the same partial differential equation [Eq.Q] with the same boundary condition [Eq.([5])]. Apart 
from their intriguing physical properties, the interest for <5-biased JTJs stems from the fact that they were successfully 
used to detect trapped fluxoids in a recent experiment aimed to study the spontaneous defect production during the 
fast quenching of a superconducting loop[l . The main peculiarity of a <5-biascd junction is the jump of the phase 
derivative at the point were the current is injected. We have shown that this discontinuity naturally leads to the 
formation of fractional vortices. In the last few years there has been a great deal of interest in the phase discontinuity 
observed in — n transition junctions which can be modeled with a bias made of two closely spaced (5-functions with 
opposite sign, more precisely by the derivative of a 5- function [3T]. We believe that our findings on the properties of 
a single <5-function can shed some more light on the mechanisms responsible of the appearance of fractional vortices 
in — 7r transition Josephson tunnel junctions. 

The authors thank A. Gordeeva for useful discussions and for helping us with the dynamical numerical simulations. 
One of us (VPK) acknowledges the financial support from the Russian Foundation for Basic Research under the grants 
09-02-00246 and 09-02-12172. 




FIG. 8: log-log plot of the measured zero-field critical currents i c (0) of many 5-biased JTJs versus their reduced length L/Xj. 
The critical current has been normalized to the small junction theoretical value calculated as the 70% of the current jump at 
the junction gap voltage. The solid line is obtained from Eq.(|ll[) with 7 C as in the empirical Eq.( |12[ ). 
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